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INTRODUCTION 

The scarcity of fresh water, the depletion of fossil fu- 
els, and the ever-increasing demand for electric power are 
important issues that receive increasing attention in a va- 
riety of branches of science and technology. In all three 
cases nanoporous carbon electrodes, immersed in a liq- 
uid medium with charge carriers, are being considered as 
device elements. For instance, in carbide-derived carbon 
electrodes with nanometer-sized pores filled with an ionic 
liquid, electric energy can be stored/released through the 
adsorption/ desorption of ionic charges on/from the sur- 
face of the charging/discharging electrodes [3-0|. Carbon 
electrodes are also being explored in capacitive devices to 
harvest sustainable energy from mixing fresh river water 
with salty sea water [4j-|7| . This salinity-gradient energy, 
or "blue" energy, is obtained from charging up a pair of 
electrodes immersed in sea water (whereby ions adsorb 
onto the electrodes at a low potential) and discharging 
them again immersed in fresh water (whereby ions des- 
orb from the electrodes at a higher potential). This ca- 
pacitive mixing process, with brackish water as a waste 
product, intercepts the spontaneous diffusion of ions from 
high to low salinity in much the same way as heat engines 
intercept the heat flow from hot to cold heat baths; for 
typical salt concentrations in river- and sea water these 
"blue engines" can produce of the order of 2kJ of energy 
per liter of river water, in principle even completely re- 
versibly Q. The reverse process, which can be seen as 
a "blue fridge", is a desalination process in which two 
volumes of initially brackish water are converted into a 
volume of fresh water and a volume of brine by charg- 
ing up the electrodes in one of the volumes (which then 
desalinates due to ion adsorption onto the electrodes, at 
a high potential) and discharging them in the other vol- 
ume (which then becomes more salty due to the release 
of the ions from the electrodes, at a low potential) 
Of course the "blue fridge" requires a net energy input, 
and ongoing research questions involve the efficiency and 
speed of such processes. 

In this contribution we will perform a thermody- 
namic and statistical-mechanical analysis of supercapaci- 
tors and blue engines. We will identify direct similarities 
between the electric work performed by/onto these de- 
vices with mechanical work performed by heat engines or 
consumed by fridges. Moreover, we will identify a num- 
ber of Maxwell relations. A distinction emerges between 



the differential capacity at constant ion number and at 
constant ion chemical potential, directly equivalent to the 
heat capacity at constant volume and constant tempera- 
ture. Finally, we will discuss the charge distribution on 
a (porous) electrode at a given potential. Throughout 
we make connection with recently published results, al- 
though we start off with a quick reminder of ordinary 
thermodynamics to clarify the analogies. 

THERMODYNAMICS OF HEAT ENGINES: A 
REMINDER 

We consider a system of internal energy U, volume 
V, and entropy S. If the number of particles in the 
system and all other (geometric, dielectric, magnetic, 
etc.) characteristics are considered fixed, we can write 
U = U(S, V) such that 

dU = TdS-pdV 1 (1) 

with temperature T and pressure p of the system defined 
by 

'-(SX-'~(3).- < 2 > 

Eq.([I]) is a combined formulation of the First and Sec- 
ond Law of Thermodynamics, where TdS is the amount 
of heat that the system takes up reversibly from a heat 
bath (also at temperature T) and pdV the reversible me- 
chanical work done by the system (on the environment 
also at pressure p). 

Consider the system to reversibly go through a cycle in 
the p—V plane, such that the final state is identical to the 
initial state. From the fact that U is a state function, we 
conclude that § dU = 0, such that the total mechanical 
work performed by the system during the cycle can be 
written as W m = § pdV = § TdS. Geometrically this 
means that the work equals the enclosed area in the p — V 
plane, but also the enclose area of the cycle in the T — S 
plane. In other words, reversible work of a cyclic heat 
engine must be accompanied by heat exchange. For the 
system to perform a positive amount of work, i.e. for the 
system to act as a heat engine, it should typically expand 
at high pressures (and hence at at high temperatures), 
thereby taking up heat from baths during (a part of) the 
expansion, and compress at lower pressures (and hence 
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at lower temperatures), thereby releasing heat into the 
colder baths during (a part of) the compression. 

Two famous examples of heat engines are the Stirling 
engine and the Carnot engine, for which the working sub- 
stance is a classical ideal gas (for which U oc T oc pV) 
that cycles in a four-fold fashion. In the (idealised) 
Stirling engine the high-temperature expansion and the 
low-temperature compression are performed isothermally 
(such that dU = and hence TdS = pdV), while 
the cooling and heating parts take place isochorically 
(dV = such that dU = TdS). In the Carnot cycle, the 
expansion and the compression consist both of an isother- 
mal part {dU = 0) and an adiabatic part (dS = such 
that dU = —pdV thereby cooling and heating the gas 
upon expanding and compressing the gas, respectively). 
It is well known that the Carnot engine yields the most 
efficient conversion of heat into work for given hot and 
cold heat reservoirs at high and low temperatures 7\ and 
Ti, respectively, with the mechanical work of the Carnot 
cycle given by W m = AT AS where AT = T h - T t and 
AS is the entropy extracted from the hot bath during 
the isothermal expansion and delivered to the cold bath 
during the isothermal compression. In the T — S plane 
this Carnot cycle is represented by a rectangular shape 
at two fixed temperatures and two fixed entropies. 



THERMODYNAMICS OF 
ELECTRODE-ELECTROLYTE SYSTEMS 

The system of our actual interest here is a macroscopic 
electrode with total charge Q in contact with a 1:1 elec- 
trolyte that contains N+Q/e counterions and N monova- 
lent coions, such that the combined electrode-electrolyte 
system is charge neutral. Here e is the proton charge. 
The temperature T is fixed, and the geometric proper- 
ties of the electrode (e.g. its surface area, its porosity, 
the volume and the curvature of its pores, etc.) are as- 
sumed to be fixed as well. Moreover, in the case of an 
aqueous electrolyte the water is treated as a structure- 
less dielectric continuum. The Helmholtz free energy of 
this system can then be written as F(N,Q), where we 
drop the dependence on the fixed temperature T and the 
fixed geometric variables for notational convenience. Re- 
gardless the functional form of F, which depends on the 
microscopic details of the ion-electrode and ion-ion inter- 
actions, we can write the differential of the state function 
F generally as 

dF = fjdN + VdQ, (3) 

where we define the ionic chemical potential p and the 
electrostatic potential of the electrode \£ as 

p=(^) and*=« . (4) 
\9Nj Q \dQ) N y ' 



We note that \i and "J, which are the intensive conjugate 
variables of the extensive variables N and Q, respectively, 
can be seen as equations of state, i.e. p = p(N,Q) and 
$ = ty(N,Q), which we leave unspecified for now as we 
focus on general and universal thermodynamic properties 
of the electrode-electrolyte system of interest. 

We first note that 9dQ in Eq.([3]) represents the 
(isothermal and reversible) electrostatic work done on the 
electrode-electrolyte system by its environment (which is 
also at electric potential 4") when it provides the sys- 
tem with an additional charge dQ. Likewise pdN is 
the (reversible and isothermal) chemical work done on 
the system when its environment (at chemical potential 
p) provides dN pairs of salt ions. Consider now a re- 
versible and isothermal cyclic process due to charging 
and discharging processes, possibly at various N, under 
the constraint that the initial and the final state of the 
electrode-electrolyte system are the same. The total elec- 
tric work performed by the system during the cycle equals 
W = — § ^dQ, i.e. the work performed is the (nega- 
tive of the) enclosed area in the — Q plane. Typically 
work is done (W > 0) by this cyclic "blue engine" if the 
electrode is charged at low voltage (which usually im- 
plies a high concentration of ions to screen the electrode 
charge) and discharged at high voltage (low salt concen- 
trations). By virtue of F being a state function, such 
that § dF = and hence W = § pdN, we indeed find 
that reversible cyclic blue engines must be accompanied 
by ion exchange processes, whereby W > implies by 
thermodynamic necessity that ions are to be taken up 
(dN > 0) by the system at a high chemical potential 
(during (a part of) the charging process) and released 
again (dN < 0) at a low chemical potential (during (a 
part of) the discharging). 

If one compares Eqs.(JT]) and © a striking resem- 
blance appears, not only regarding the (free) energy con- 
tribution due to the exchange of heat and ions, TdS 
and pdN, but also regarding the mechanical and elec- 
tric work contributions, —pdV and ^dQ, respectively. 
In fact, we can make the following mapping: U o F, 
T O p, S o N, -p o and V o Q. This iden- 
tification of variables respects the symmetry of exten- 
sivity and intensivity, connects (de)compressions of a 
gas with (dis) charging of the electrolyte-immersed elec- 
trode, and implies that isothermal (constant-T) and adi- 
abatic (constant-.!?) volume changes in heat engines are 
analogous to grand-canonical (constant-/^) and canonical 
(constant- N) charging processes of electrolyte-immersed 
electrodes, respectively. 

With this mapping of variables between heat engines 
and blue engines in mind, it is interesting to note that 
the blue engine recently developed by Brogioli [3| to har- 
vest salinity gradient energy is actually equivalent to 
the (idealised) Stirling engine: Brogioli's electrode charg- 
ing/discharging processes take place at constant p, just 
like Stirling's volume changes take place at constant T, 



3 



and Brogioli's flushing processes to exchange river- and 
sea water (to change the ion chemical potential) at con- 
stant electrode charge are equivalent to Stirling's heat 
exchanges between hot and cold baths (to change the 
temperature) at constant volume. In addition, this map- 
ping of variables was exploited in Ref.[|| to construct a 
conceptual Carnot-like blue-engine. The key difference 
with the Brogioli cycle is the replacement of the two 
flushing steps by a constant-iV charging and discharg- 
ing process, whereby the initially salty water desalinates 
upon electrode charging due to ion adsorption, and the 
initially fresh water salinates upon discharging due to ion 
desorption. This Carnot-like cycle composed of two iso-/z 
and two iso-N (dis)charging steps, yields an electric work 
output W — ANAfj, per cycle, where A7V is the number 
of ion pairs that flows during the iso— /i parts of the cycle 
from the salty to the fresh water through an adsorption- 
desorption process onto the electrodes, and where A/i is 
the chemical potential difference between the salty and 
the fresh water. The rectangular shape of the enclosed 
area of the cycle in the /j,-N representation is the hall- 
mark for the most efficient process, since it makes explicit 
that each of the AN ion pairs contributes its full chemi- 
cal potential difference Afi to the total work; more work 
(per cycle per transferred ion pair) is thermodynamically 
impossible. 

There is, however, one key difference between the na- 
ture of the mapped variables T and fj,, since T does have 
a well-defined absolute zero whereas /i is only defined up 
to an arbitrary reference potential. For that reason there 
is no well-defined analogue of the Carnot heat-engine effi- 
ciency AT/Th', the analogous expression Afi/fih with jjih 
the high chemical potential of the sea water is meaning- 
less. 



MAXWELL RELATIONS AND RESPONSE 
FUNCTIONS 

Starting from Eq.Q a number of of thermodynamic 
relations for (super)capacitors immersed in an ionic fluid 
can be constructed in full analogy to the standard re- 
lations that follow from Eq.([T]) for heat exchange and 
volume work. In the latter case it proves convenient, 
for instance, to consider Legendre transformations of 
U(S, V) to obtain thermodynamic potentials such as the 
Helmholtz free energy, the Gibbs free energy, or the grand 
potential, from which a set of Maxwell relations follows 
by considering second derivatives with respect to two 
different variables What emerges from this is a 

set of five measurable and therefore interesting response 
functions associated with heat transfer and pressure- 
volume work: the constant- volume and constant pressure 
heat capacities c v = T{dS/dT) v and c p = T{dS/dT) p , 
the isothermal and the adiabatic compressibilities kt = 
-V- x (dV/dp) T and k s = -V- 1 {dV/dp) Sl and the iso- 



baric thermal expansivity j3 p = V~ 1 (dV/dT) p . However, 
standard textbook thermodynamics dictates (on the ba- 
sis of reciprocity relations and Maxwell relations) that 
these five quantities are not all independent, as they sat- 
isfy 
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and 



cv 



T 



(5) 



such that there are in fact only three independent re- 
sponse functions. As a consequence of the second rela- 
tion together with n T > 0, it is guaranteed that c p > c v 
and hence from the first one that > Kg. Below we fol- 
low exactly the same thermodynamic arguments for the 
electrode-electrolyte system of interest here, building on 
the mapping of the variables as discussed in the previous 
section. 

We consider three Legendre transformations of 
F(N,Q), which we denote by the Gibbs-like free en- 
ergy G(N, = F - *<3, the grand-potential ft(fi, Q) = 
F — fiN, and the thermodynamic potential Y(n, = 
F — fiN — *ifQ. Note that all these potentials depend im- 
plicitly also on T and on the geometric properties of the 
electrode (e.g. the volume), which renders even Y well- 
defined as any bonafide thermodynamic potential must 
depend on at least one extensive variable. The differen- 
tials of these potentials are given by 



dG = fidN-QdV, 
dQ = -Ndfi + ^dQ, 
dY = -Ndfi-QdV. 



(6) 



By taking "off-diagonal" second derivatives of each of the 
four potentials, the following four Maxwell equations can 
straightforwardly be derived: 



d 2 F 
dNdQ 

d 2 G 
dNd^ 

d 2 n 

dfidQ 
d 2 Y 



dN 



9Q/ N 
dp, 

( dN 
dQ\ 



(7) 
(8) 



We note that an alternative derivation of Eq.© was 
recently reported, and in fact both sides of the equa- 
tion as obtained from measurements were successfully 
compared 0, where the explicit ideal-solution relation 
between chemical potential and salt concentration was 
used. We also note that as defined in Eq.© plays 
the same role here as the isobaric expansivity j3 p in "stan- 
dard" thermodynamics. 

By considering "diagonal" second derivatives of the 
two potentials that depend on Q we define the iso-/z and 



4 



the iso-A capacitances Cn and Cn 



a 2 Fy 1 


(dQ 


dQ 2 ) 






(dQ 


dQ 2 ) 





N 



(9) 
(10) 



which we recognise on the basis of our mapping of vari- 
ables as the analogues of the two compressibilities. Tak- 
ing "diagonal" second derivatives of the the two poten- 
tials that depend on N we find 



d 2 F\ _1 


(dN 


dN 2 ) 
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(dN 


dN 2 ) 





XQ 



X*, 



(11) 



where the x's play the role of the heat capacities. 

The five quantities C^ , Cn, x* ; XQ; and are 
not all independent, and satisfy relations akin to the two 
standard relations between c p , cy, Kg, kt: and P p of 
Eq.©. Standard thermodynamic manipulations involv- 
ing the reciprocal and the reciprocity relations yield 



Cn 



XQ 



— and 



X* 



a* 

XQ = 7T- 



(12) 



With C M > 0, which is a stability requirement as we will 
see below, we thus find that > XQ and C M > Cn- 
Note that thermodynamics does not provide numerical 
values for these quantities, as this would require a micro- 
scopic or molecular theory for the electrode-electrolyte 
system. The strength of these thermodynamic relations 
lies, however, in their generality: it is thermodynami- 
cally guaranteed for any electrode-electrolyte system that 
C^ > Cn- It is therefore guaranteed that *(Q) rises 
faster with Q at fixed N than at fixed /j,, the difference 
being larger if a* is larger, i.e. if the electrode exhibits a 
larger adsorption-to-potential or charge-to-concentration 
response. Given that a large difference C M — Cn gives 
rise to a large area in the enclosed \t — Q plane, and 
hence a large amount of work during a (Carnot-like) cy- 
cle, it could be beneficial for these devices to be based on 
electrode-electrolyte combinations with a large a^. 



ENSEMBLES AND CHARGE DISTRIBUTION 

We have been concerned with an electrode-electrolyte 
system of which the thermodynamics can be described 
by the Hclmholtz free energy F(N, Q). This system has 
an underlying microscopic Hamiltonian, denoted by H 
here, that depends on all the degrees of freedom of all the 
ions and all the surface charges. At temperature T, the 
statistical probability of a microscopic configuration is 
then given by the Boltzmann weight exp(—/3H)/Z(N, Q), 



where the normalization factor Z(N, Q) is the canonical 
partition function and where /3~ 1 = k B T . The Helmholtz 
free energy, with the differential given in Eq. (J3j, follows 
as F(N,Q) = -k B TlnZ(N,Q) as usual. 

If we now consider the electrode-electrolyte system at 
fixed N and the total electrode charge Q is a fluctuat- 
ing quantity that takes values according to the thermal 
probability distribution 



P(Q) = 



exp(-/3F(JV, Q) + 
Z(N, *) 



(13) 



where the normalisation factor is the Gibbs-like parti- 
tion function Z(N, *) = £ Q exp(-/9F(iV, Q) + /3#Q) = 
exp(-/3G(iV, *)), with G{N,V) the Gibbs-like potential 
with a differential given by Eq. (J6]) . It readily follows 
that the average electrode charge is given by {Q)n = 
Eq P (Q)Q = k B TZ- l {dZ/d^) N = -(8G/d9) N , in 
agreement with the differential of Eq.©. Likewise one 
can write (Q 2 ) N = {k B Tf Z~ l (d 2 Z / d^ 2 ) N , which af- 
ter some elementary algebra gives rise to the following 
identity for the variance 



(Q 2 )n - (Q) N = k B TC 



(14) 



with Cn the constant- N differential capacity defined in 
Eq.©. One can also show that the variance of Q at 
fixed fi and * is given by k B TC^, which is larger than 
in the iso-7V case since C^ > Cn as we have seen above. 
Expression (fT4"l) shows that the differential capacitance 
can be measured from the fluctuations of the electrode 
charge, completely equivalently to measurements of the 
heat capacity from energy fluctuations and the compress- 
ibility from volume fluctuations in "ordinary" (NVU) 
and (NVT) ensembles, respectively. 

Of course the differential capacity C (either Cn or 
C M ) as well as the average electrode charge (Q) (either 
(Q)n or (Q)fj,) are extensive quantities that scale lin- 
early with the system size (in this case the electrode 
area). Therefore the standard deviation of the charge, 
(ksTC) 1 / 2 = SQ, becomes much smaller than the aver- 
age charge (Q) for thermodynamically large electrodes. 
However, our analysis can also be applied to a small 
(sub-)system provided it is large enough to be statisti- 
cally independent (in fact F(N, Q) = F(NM,QM)/M 
for M > 1 must hold). An example includes typical 
computer simulations of 10 2 -10 4 ions near electrode ar- 
eas of the order of tens of nm 2 . For such small electrode 
areas the charge fluctuations at fixed \& can be signifi- 
cant, and in fact even be of the order of (Q). Consider, 
for instance, a typical carbon-based supercapacitor with 
an areal capacity of the order of several /iF/cm 2 at a 
potential W = IV, such that the average charge-density 
is of the order 10 _2 e/nm 2 . A patch of electrode of the 
order of 100 nm 2 contains, therefore, a charge (Q) ± SQ 
of the order of e ± e, indicating that a significant fraction 
of nm-sized patches carries a charge that is opposite to 
the average charge. 



5 



In fact it is possible to calculate the complete charge 
distribution P(Q) as defined in Eq. ([T3f from the formal- 
ism that we used here by expanding F(N, Q) about the 
most-probable charge Q* , defined by P'(Q*) = 0, i.e. by 
F'{N, Q*) = Vf, where a prime denotes a derivative with 
respect to Q. Skipping the iV-depence for notational con- 
venience, we then find 

P(Q) cx exp(-i(Q-Q*) 2 F"(Q*) 

( (Q - Q*f (Q-Q*fC'(Q*) \ 
6XP v 2k B TC* 6k B T(C*) 2 ' ' ' J ' 

(15) 

where C* = Cn(Q*)- On the basis of extensivity argu- 
ments the cubic and higher-order terms in the exponent 
can be ignored in the thermodynamic limit, such that 
P(Q) is a Gaussian with Q* — (Q) and a variance in ac- 
cordance with Eq. (fT4")) . However, for smaller systems the 
higher order terms may be relevant, at least in the case 
when C'(Q) = (dC/dQ)^ 7^ 0, i.e. when the differen- 
tial capacity depends significantly on the average charge 
(and hence on the applied voltage) [Tljj . Eq.(fl~5"|) appears 
to be in agreement with the findings in Ref. 3] on graphite 
electrodes, which give an almost Gaussian charge distri- 
bution and an essentially vanishing C. However, the 
significantly skewed charge distributions as found Ref. 3] 
for nanoporous carbide-derived carbon electrodes require 
some further attention, as the capacity is reported to be 
essentially constant, which should yield a Gaussian distri- 
bution according to the present derivation. We speculate 
that the individual carbon atoms of the electrode may 
be too small to be viewed as a statistically independent 
subsystem, although more research is needed to clarify 
this point. 

SUMMARY 

We compare the differential of the Helmholtz free en- 
ergy F(N, Q) of an electrode-electrolyte system with that 
of the energy U(S, V) of an "ordinary" thermal system, 



and identified the transferred heat TdS with the ion 
flow contribution /j,dN, and the mechanical work —pdV 
with the electric work *S>dQ. By a mapping of the vari- 
ables (N, Q,^) of present interest onto (S, T, V, — p), 
we can define Legendre transformations of F(N, Q), iden- 
tify Maxwell relations, and formulate Eq. dT^I) analogous 
to ([5]). With this mapping we identify the Brogioli blue 
engine as a Stirling heat engine, and we discuss a Carnot- 
like blue engine that should have the optimal conversion 
of mixing entropy to work. Finally we discuss the charge 
distribution at fixed electrode potential, and show that 
the variance of the charge scales with the differential ca- 
pacitance. 
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